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Abstract 


The RICE experiment (Radio Ice Cherenkov Experiment) at the South Pole, 
co-deployed with the AMANDA experiment, seeks to detect ultra-high en¬ 
ergy (UHE) electron neutrinos interacting in cold polar ice. Such interactions 
produce electromagnetic showers, which emit radio-frequency Cherenkov ra¬ 
diation. We describe the experimental apparatus and the procedures used to 
measure the neutrino flux. 
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I. ULTRA-HIGH ENERGY NEUTRINO ASTRONOMY: INTRODUCTION 


Detection of ultra-high energy [Ey > lO^^eV) neutrinos represents a unique opportunity 
to probe the distant universe. High-energy protons and photons from distant sources are 
likely to interact with the cosmic microwave background (CMB); protons, being charged, 
have their trajectories bent in galactic and intergalactic magnetic fields. Neutrinos are inert 
to CMB photons and point directly back to their source, giving essential information on 
those sources. In the realm of particle physics, detection of UHE neutrinos from cosmologi¬ 
cal distances, if accompanied by flavor identification, may permit measurement of neutrino 
oscillation parameters over an unprecedented range of Am^. It has been suggested that, 
with a sensitive enough array, tan neutrinos may be identified by two-step “double-bang” 
|jl[] processes where a r lepton is created and subsequently decays. Since neutrino absorp¬ 
tion in the earth depends on the chord length through the earth,Q the angular distribution 
of detected neutrino events could be used to verify predictions for weak cross-sections at 
energies unattainable by any man-made accelerator. Alternately, if the high energy weak 
cross-sections are known, they can be used to test Earth composition models along an arbi¬ 
trary chord (‘neutrino tomography’) |^. 


A. Current Experimental Efforts 

Several recent projects (including AMANDA 0, NESTOR 0, Lake Baikal [0, 
ANTARES PI) are optimized for detection of very high energy (10^^“^® eV) cosmic ray 
muon neutrinos. Sensitivities to higher energies, as well as electromagnetic cascades, have 
also been shown to be substantial in such experiments p,|[. These instruments are based on 
photomultiplier tube detection of the optical Cherenkov cone from muons produced in muon 
neutrino charged current interactions. At high energies, muons have ranges of order 1 km 
and follow approximately straight trajectories (smeared by multiple scattering), punctuated 
by catastrophic bremstrahlung every 0.1-1 km or so, in which ~10% of the muon’s energy 
is lost to a photon. RICE employs radio detection, which is believed to be the most efficient 
detection mechanism at energies of 10^® eV and beyond p. 

The RICE concept is illustrated in Figure |^, depicting a Cherenkov cone ht to a set of 
“struck” dipole receivers in a simulation event, along with the extracted neutrino direction. 
Receiver locations are drawn to scale in the Figure. In the actual array geometry, dipole 
receivers are spread over a 200 m x 200 m x 200 m cube beneath and around the Martin 
A. Pomerantz Observatory (MAPO), approximately 1 km. from the geographic South Pole. 


^The Earth is opaque to >PeV neutrinos with Standard Model cross-sections and with zenith 
angles approaching 180 degrees. Largely because of earth shadowing, the RICE array is most 
sensitive to neutrinos incident at zenith angles between 60 and 120 degrees. Conversely, the angnlar 
distribution is sensitive to the cross-section and may allow checks of the Standard Model. 
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FIG. 1. Simulated RICE event. The actual detector geometry is shown, to scale. 


B. Radio Detection 


We have initiated a pilot program based on radiowave receiver technology in order to 
extend electron neutrino detection up to the PeV energy scale. The detector is intended 
to have sensitivity in the energy range ~ 10^® — 10^^ eV. RICE (“Radio Ice Cherenkov 
Experiment”) should ultimately be capable of observing the sky with angular resolution of 
10-50 mrad. 

Coherent radio Cherenkov emission is an efficient method for detecting high energy par¬ 
ticles. The history of the effect goes back to Jelley, who hrst considered whether cosmic 
ray air showers might produce a radio signal Askaryan [^] subsequently predicted a 


net charge imbalance in air showers, and coherent radio power proportional to the energy of 
the shower squared. Substantial radio emission from atmospheric electromagnetic cascades 
was observed more than 30 years ago Progress in ultra-high energy air showers has 


sparked renewed interest, and new observations of radio pulses have been reported recently 
m , suggesting a possible radio component to the Auger detector [jT^. A recent international 


meeting highlights current progress 

This effect has recently been observed in a test-beam experiment at SLAC [0 . A beam 
of electrons, of known amperage, was hred into a sand target and the radiation resulting 
from the impact measured in the radio regime. As Figure illustrates, the measured signal 
strength (diamonds) displays the expected dependence on both total beam current (left, 
solid curve) and frequency (right, dashed curve). 


C. Initiation of the RICE Experiment 

The RICE experiment was initiated in October, 1995 when the AMANDA collaboration 
graciously consented to co-deployment of two shallow radio receivers (“Rx”) in the hrst holes 
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FIG. 2. Results of testbeam experiment[14], showing signal strength and expected dependence 
on total beam current (circles, left) and frequency (diamonds, right). Figure reprinted courtesy of 
D. Saltzberg and P. Gorham. 

being drilled for AMANDA-B. Following deployment, a surface transmitter (“Tx”) was used 
to verify that signals could be detected under-ice with better than 10 ns timing precision. 
However, cross-talk and ampliher oscillation problems precluded use of those receivers for 
science. The hrst three dedicated RICE receivers were deployed in 1996-97, along with one 
underice transmitter. The extraordinary RF clarity of South Polar ice was amply illustrated 
by the evident brightness of the AMANDA photomultiplier tubes 2 km. below the RICE 
receivers (and close to the null in the dipole antenna’s reception pattern). A Fourier analysis 
of the RF transients produced by AMANDA photo-tubes in laboratory conditions indicated 
that PMT background power dominated at frequencies below 100 MHz, motivating use of 
high-pass hlters in subsequent Rx deployments. Antenna deployments followed in 1997-98 
(three more receivers and two more transmitters deployed in AMANDA-B holes), 1998-99 
(six receivers deployed in 5”-diameter dedicated RICE holes, bored using a mechanical, 
rather than a hot-water drill), and 1999-2000 (six receivers, and one transmitter deployed in 
AMANDA-B holes). 


II. EXPERIMENTAL APPARATUS 

The RICE experiment presently consists of an 18-channel array of radio receivers (“Rx”), 
scattered within a 200 mx200 mx200 m cube, at 100-300 m depths. Each receiver contains a 
half-wave dipole antenna, offering good reception over the range 0.2-1 GHz. Twelve receivers 
are buried in the boreholes drilled for the AMANDA photomultiplier tube deployment during 
the 1996-97, 97-98, and 99-00 austral summers. Six receivers are located in dedicated RICE 
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holes; four such holes were drilled with a 5-inch diameter mechanical hole-borer in 1998- 
99. The signal from each antenna is immediately boosted by a 36-dB in-ice ampliher, then 
carried by ~300 m coaxial cable to the surface observatory, where the signal is hltered 
(suppressing noise below 200 MHz due to both AMANDA photo-tubes, as well as continuous 
wave backgrounds from South Pole station at 149 MHz), re-ampli£ed (either 52- or 60-dB 
gain), and split into two copies. One copy is fed into a CAMAC crate from which, after 
initial discrimination (using a LeCroy 3412 discriminator), the signal is routed into a NIM 
crate where the trigger logic resides. The other copy of the analog signal from the antennas is 
input to one channel of a digital oscilloscope, where waveform information is recorded. Also 
deployed are three large TEM surface horn antennas which are used as a veto of surface¬ 
generated noise. 


A. Detector Array Status 

The status of the current array deployment is summarized in Figure ^ Indicated in 
the Figure are the AMANDA holes (1-19, drilled using the hot-water technique) containing 
RICE receivers, and also the four holes drilled specifically for RICE in 1998-99 containing 
RICE-only equipment. All channels are indicated by a 5-character alphanumeric mnemonic 
corresponding to the year of deployment, the type of dipole antenna deployed (transmitter 
“Tx” or receiver “Rx”; these only differ by the installation of a receiver ampliher in the latter 
modules), and a numerical identiher. Also indicated in the Figure is the MAPO building 
which houses hardware for several experiments, including the RICE and AMANDA surface 
electronics. The AMANDA array is located approximately 600 m (AMANDA-A) to 1400 
m (AMANDA-B) below the RICE array in the ice; the South Pole Air Shower Experiment 
(SPASE) is located on the surface at a; ~ —400m in the Figure. The coordinate system 
conforms to the convention used by the AMANDA experiment; grid North coincides with 
the -l-y-direction in the Figure. 


B. Event Trigger and Vetoes 

The main RICE trigger condition requires that four or more channels have their output 
voltage exceed a common adjustable threshold within a pre-set time window (currently 
At = 1.2/is, set by the radio wave transit time across the array). The threshold is adjusted 
so that thermal huctuations, or other backgrounds, do not cause excessive triggers. The 
time window allows for coincidences across the full array, independent of event geometry. 
A valid event trigger is also dehned when, in a time window of 1.2 ps, i) > one underice 
antenna registers a signal above threshold in coincidence with a 30-fold PMT AMANDA-B 
trigger, or: ii) > one underice antenna registers a signal above threshold in coincidence with 
a high-multiplicity SPASE event. 

Additionally, there are two ways that surface-generated background transient events can 
be vetoed - either: a) one of the surface horn antennas registers a signal above threshold, in 
which case data-taking is inhibited for the subsequent 4 /isec,0 or b) the timing sequence of 


^The raw rate for this veto is of order 10 Hz, so this represents an insignificant loss of data. 
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FIG. 3. Present geometry of the RICE array, relative to AMANDA hole 4. “Tx” designate 
transmitters; “Rx” designate receivers. Depths of receivers, as well as relative cable diameter 
(indicated by the thickness of horizontal lines before the antenna identifier), is also indicated. 


6 






























































hits in the underice antennas is determined to be consistent (in software) with the sequence 
expected from surface-generated backgrounds.^ This “surface-veto” algorithm uses the time 
differences between TDC hits recorded on-line to make a fast (<10 msec) classification of 
the event as surface/non-surface in origin.^ One of every 10000 events classified as “surface” 
(i.e., having a sequence of antenna hits consistent with a Z >0 vertex) are retained to ensure 
that this veto is functioning properly. 

If any of the above trigger criteria are satisfied and there is no veto signal, the time of 
each hit above threshold (as recorded by a LeCroy 3377 TDC), and also an 8.192 ps buffer 
of data stored in an HP54542 digital oscilloscope at 1 GSa/s (for each channel) is written to 
disk. Each event is also given a GPS time stamp for synchronization with other South Pole 
(and more global) experiments. To monitor the stability of the amplifiers as well as changing 
background conditions, waveform measurements are taken every ten minutes, independent 
of other event triggers (so-called “unbiased” triggers). A schematic of the event trigger is 
shown in Figure 


1. Discriminator Efficiency 


As the first element in the event trigger, it is essential that the Lecroy 3140 Discriminator 
module operates efficiently for nanosecond duration pulses. The discriminator efficiency is 
checked explicitly in the lab using an 800 ps width signal (generated by an HP8133A signal 
generator) fed directly into a 3140 Discriminator module, and also verified in the field by 
pulsing the array with a 600 ps width signal broadcast through one of the broadband TEM 
horns (the bandwidth of the horns is ~3 GHz, so timing is preserved down to ~300 ps). 
The LeCroy discriminator is found to be >99% efficient in triggering on the narrowest width 
pulses our signal generator is capable of producing. Given the measured impedances of our 
cable, amplifiers, and antennas (discussed below), a neutrino pulse is expected to produce 
signals of width ~l-3 ns at the surface. 


C. Raw Data 

Raw trigger rates (before veto) are typically 30 Hz; data-taking rates after the veto 
are typically 0.01 - 0.1 Hz. The data-taking rate (and, correspondingly, our experimental 
livetime) is limited by: a) the time required to write information from the digital oscilloscopes 
to disk (~10 sec/event), b) the time required to perform the surface-background veto in 
software (~10 msec/event), and c) our inability to take data at those times when the South 
Pole satellite uplink (broadcast at 303 MHz) is active, due to the high amplitude of the uplink 
signal.^] When this uplink is not active, the discriminator thresholds correspond to typical 


^In the event of an AMANDA or SPASE coincidence trigger, the surface veto is disabled in order 
to preserve possible air shower events. 

^This functionality will reside in a dedicated CAMAC board beginning in 2002. 

^Notch filters will be installed in 2002 to mitigate this background. 
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FIG. 4. Functional schematic of RICE event trigger. 




















































































livetimes of ~80%. Based on known dead times in the system, integrated array on-time is 
cnrrently monitored and npdated on-line. 

A typical event, as recorded in the HP54542A oscilloscopes, is shown in Fignre In the 


edge-c-big50-6-2-11-99_10.dat 988.76083 Veto 



FIG. 5. Typical data event recorded by the RICE array (Feb. 11, 1999). This event was vetoed 
as having a surface origin, based on the observed sequence of hit times (see text). 

event display, each of the 8.912 /rsec traces corresponds to one of the under-ice antennas. 
Vertical bars indicate the TDC time recorded for each particular channel. Note that one 
channel (12) has no waveform information, and two channels are not used in the trigger 
(13 and 15) and therefore have no TDC information. Also note that channels which do not 
have hits which exceed the initial discriminator threshold do not register a TDC time. The 
rise time and leading-edge timing resolution of each receiver as determined by waveform 
information can be approximately estimated from the Figure. The leading edge resolution 
is typically 2 nsec; the ring time for each antenna is typically tens of nanoseconds. In the 
Figure, the antennas have been ordered according to their distance from the surface, with 
the shallowest (/deepest) antennas at the top (/bottom) of the Figure. The pattern of hits 
shown is therefore easily identihable as a surface-generated RF pulse sweeping down through 
the radio array. 

The Fourier transform of a typical waveform signal (channel 7) is displayed in Figure ^ 
The general shape of the signal in the frequency domain reveals many interesting features of 
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the experimental hardware and the environment. Below 200 MHz, suppression of noise due 
to the filter is evident; as mentioned previously, use of this hlter was motivated by the need to 
reduce low-frequency RF backgrounds due to broadband galactic sky noise, monochromatic 
Continuous Wave (CW) sources at the Pole, as well as the firing of AMANDA phototubes. 
Above 200 MHz, losses due to cable effects become increasingly pronounced. At the time 
this event was recorded, the 303 MHz satellite uplink at the Pole was active;^ the peak at 
this frequency is evident in this Fourier Transform. 
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FIG. 6. Fourier transform of raw data waveform corresponding to channel 7 of the previous 
Figure. The South Pole station satellite uplink at 303 MHz is evident in the Figure. A high pass 
filter suppresses signals below ~180 MHz. Other features of this distribution are discussed in the 
text. No corrections for cable loss, signal amplification, etc. have been made. (The vertical scale 
is arbitrary.) 


Spectrum for channel 7 98Rx2 



III. TIMING CALIBRATION 

Event and source reconstruction is based on our knowledge of the array geometry, ice 
properties and thus the expected times for a wave-front to propagate from the source to any 
given receiver location. A time-of-hit is defined as the time of the hrst excursion exceed¬ 
ing 6arms in a waveform; arms is determined from the sample of “unbiased” events which 
are taken every 600 seconds, independent of the event trigger status. The maximum reso¬ 
lution on the hit-times cannot, therefore, exceed the sampling time of the oscilloscopes (1 


®This event was taken before our “303-veto” was functional. 
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ns). This resolution can be improved using a “matched hlter” algorithm which matches 
the observed waveform with a reference signal waveform; such a software algorithm is cur¬ 
rently under development. In the “matched hlter” approach, one dehnes a reference signal 
corresponding to the expected damped oscillator response of the antenna to an impulsive 
signal: f{t) ~ cos {out), with a; = RC. The values of R, L and C can be estimated 

from network analyzer measurements of the complex antenna impedance Zantenna (described 
below), as a function of frequency. The signal time is dehned as that time which maximizes 
the integral: / f{t)V{t)dt, with V{t) given by the direct data measurement of the waveform 
voltage with time. As an illustration of the potential improvement arising from the matched 
hlter algorithm. Figure compares the raw signal V{t) (top panels) with the output of the 
hlter algorithm (lower panels) for one channel (channel 2) in one event. The signal to noise 
is clearly superior in the latter case. 



FIG. 7. Raw V(t) for two waveforms (top); the signal time as determined by a simple maximum 
voltage algorithm corresponds to the maximum excursion from zero. Bottom displays show the 
quantity T,f{t)V{t), with f{t) being the damped oscillator reference signal and V{t) the actual 
recorded waveform shown in top panels. The signal to noise is clearly greater for the filter output. 

Knowing the time diherences 6tij between all pairs {i,j) of hit antennas, we perform 
minimization to hnd a source location. Once the vertex location has been found, superim- 
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posing a Cherenkov cone of half-width 57°[| on the hit antennas allows a determination of the 
source direction. For the full reconstruction, this method requires at least four antennas to 
be hit (i.e., SSUj values). Uncertainties in 6tij arise from several sources, including risetime 
resolutions (~2 ns), differences in signal propagation velocity in the ice due to variations in 
the dielectric constant with depth, differences in signal propagation speed within the differ¬ 
ent analog cables being used, differences in cable lengths, and receiver deployment surveying 
uncertainties. 

Buried transmitters (“Tx”) are used to calibrate the channel-to-channel timing delays. 
A short duration pulse is sent to one of the five under-ice transmitters, which subsequently 
broadcasts the signal to the receiver array. An event vertex is reconstructed exclusively from 
the measured channel-to-channel time delays; constraining the source to a unique location 
allows a calculation of the timing residual each channel, based on the timing un- 

^^measured_^^expected 

certainty: xliij) = iterative procedure is used to calibrate out the 

observed channel-to-channel timing delays and minimize the timing residuals for an ensemble 
of events. Typical timing calibration corrections are ~10 ns per channel; these corrections 
are then used for all subsequent event reconstruction. 

Once the source location has been determined, the timing resolution for each channel 
can be derived by examining the width of the distribution. After time 

calibrations have been performed, these differences of time differences should (ideally) be 
zero for a unique vertex. Figure ^ shows these distributions, with channel 12 (arbitrarily) 
chosen as the reference channel. From these plots, the signal arrival time resolution for good 
channels (i.e., all channels except channel 8 in the Figure) is determined to be ~1.5-2 ns, 
slightly larger than the oscilloscope sampling time of 1 ns. 


A. Vertex Reconstruction Algorithms 

Two algorithms are used to determine a source location. The first, a “trial-and-error” 
procedure, tests the consistency of each possible (x,y,z) source location, in a grid 2km x 2km x 
1km below and around the radio array, with the observed data, he., the {xgrid, ygrid, zgrid) 
source point (using a step size of 10 m, or a total of 4 x 10® source points tested) corresponding 
to the minimum is identified. The advantage of this procedure is that time distortion 
effects due to ray tracing through regions with varying index of refraction can easily be 
incorporated; the disadvantage is speed and the intrinsic dependence on the grid spacing. 
The second procedure analytically determines the source location rsource and a global event 
time to by simultaneous solution of the 4 equations: \rRx,i — 'rsource\ = {,c/n)ti {i = 1, 2, 3,4), 
where is the vector from the origin to receiver i, rsource is the vector from the origin 
to the source, and ti is the hit time recorded for the receiver. There are at most two 
roots: eliminating complex or acausal roots, while also requiring consistency of 5 or more 
hits resolves ambiguities. In Cartesian coordinates, we denote this analytic solution for the 
source location as fsource={x4:hit,y4:hit, z-ihit). 

Source reconstruction results are shown in Figure Three plots, each showing the 
difference between the reconstructed vertex location and the surveyed vertex location in one 


^nice(200 MHz)~1.78. 
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FIG. 8. Time difference between calibration signal arrival time in ch. i relative to ch. 12, in 
nanoseconds. Data taken from transmitter pulser event sample. Oscilloscope jitter in channel 8 is 
evident from the distributions. 

of the spatial coordinates, are shown. Based on the (general) consistency between both the 
4-hit and the grid algorithms with the surveyed Tx location, we conclude that the surveying 
error in the Tx location is less than 5 m. 

For general events (such as those displayed in Figure we expect backgrounds generated 
at the surface to dominate. Figure 0 shows the reconstructed source depth {zgrid) for a 
sample of general triggers; as expected, sources populate the region above the array.0 


B. Dielectric Constant of Ice 

The complex dielectric constant e(a;) allows one to calculate both the absorptive and 
refractive effects of the ice. The imaginary part of the dielectric constant (related to the 
“loss tangent”) prescribes losses due to absorption; the real part corresponds to the refractive 
index of the medium. Information on the refractive index can be derived from temperature 
and density profiles, as a function of depth, acquired by AMANDA deep drilling operations. 
Such profiles {z{T,p)) can be combined with laboratory measurements of the dependence of 
the index of refraction of ice on temperature and density (n(T, p)) to predict the expected 
index of refraction profile at the South Pole as a function of depth n{z). This derived n{z) 


®Note that the z = 0 bin includes all source depths corresponding to z > 0, as well. 
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FIG. 9. Reconstruction of Tx pulse from transmitter 97Tx2, based exclusively on measured 
receiver times. 4-hit solutions are shown as the unshaded histogram; grid solutions to the vertex 
location are shown as the shaded histogram. 


function can be compared with the prohle calculated directly from RICE transmitter data. 
Using the surveyed transmitter and receiver locations, combined with Fermat’s principle, one 
can determine the profile n{z) which best reproduces the observed Tx —>■ Rx radio signal 
transit times. A preliminary comparison between the measured n{z) function determined 
from transmitter data with the “derived” n{z) profile is shown in Figure [^. Agreement is 
satisfactory; however, work is still in progress to quantify the agreement between the two 
curves. 

Absorption of radio waves in cold glacial ice is temperature, density, and frequency depen¬ 
dent and has been measured [0, indicating attenuation lengths a ~1-10 km at frequencies 
in the range 100 MHz - 1 GHz. Given the small scale of our current array (~100 m) com¬ 
pared to the very large attenuation length a expected for Polar ice in the 100 MHz - 1 GHz 
regime, our transmitter tests are consistent with zero absorption over our short baseline; 
long baseline transmitter measurements are an objective of future Polar campaigns. 
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Depth distribution for general triggers (minimal cuts) 



FIG. 10. Typical depth distribution for a sample of 348 general RICE triggers. Triggers are 
observed to be dominated by sources near the surface. (There are no underflow entries.) 
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FIG. 11. Index of refraction as a function of depth. 


IV. ANTENNA RESPONSE AND GAIN CALIBRATION 

To detect a neutrino of a given energy interacting at a given distance from an antenna, 
we need to quantify the sensitivity of the array to the resulting Cherenkov signal. Although 
neutrino vertex location reconstruction is based only on channel-to-channel timing, an am¬ 
plitude calibration is needed in order to ensure that the discriminator efficiency, for a given 
incident neutrino energy, is reliably calculable. 


A. Antenna Effective Area 

One of the fundamental parameters used to dehne the power response of an antenna 
is the effective area Aeff{0,(f)) (or gain G = iTiAeff/X"^); this is related to the incident 
signal intensity / as: Pout = We have used two techniques to measure the effec¬ 

tive area A^ff of RICE dipole antennas. The hrst applies the Friis Transformation Equa¬ 
tion to data taken when a transmitter (Tx) broadcasts to a receiver (Rx) on the KU An¬ 
tenna Testing Range (KUATR). The Friis Equation relates the power broadcast by the 
transmitter to the signal intensity measured at the Rx (valid in the far-held case). Dehn- 
ing the monochromatic power into the transmitter PTx,in, the power broadcast outwards 
from the transmitter PTx,out, the power intercepted by the receiver PRx,in, and the power 
transmitted out the back of the receiver: PRx,out (which we measure either on an oscil¬ 
loscope or a network analyzer), we have: PTx,out = Ptxm^Tx, Irxm = PTx,out/(,4:7iR‘^), 
PRx,out IRx,inAeff^Rx PTx,outAoff^Rx /)j TTeTT ^TxA^ff^Rx /). Usiug a Cali¬ 
brated transmitter, Gtx is known; by comparing to a calibrated receiver known), 

the gain or effective area of any arbitrary antenna can be inferred by a simple ratio. 

Alternately, the gain/effective area can be determined by sending a known amount of 
power through a cable and into an antenna. The rehection coefficient Tr (having real and 
imaginary components T^ and Tj, respectively) measures the impedance mismatch of the an¬ 
tenna vecZantenna and the cable Zcable, and is ideally zero for a high-gain antenna (impedance 
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well-matched to 50f2 cable, e.g.). Numerically, G ~ (1 — iF/jp). The two measurements of 
Aeff (from the testing range and the impedance measurement, respectively) agree to within 
±1 dB over the frequency range of interest and are consistent with simple expectations for 
dipoles. 


B. Effective Height 


To understand the shape of the signal voltage produced by a receiver in the time domain, 
we must quantify the complex effective height h The effective height h (in units of 


meters) is related to the magnitude and the phase of the voltage resulting from the application 
of a complex electric held vector at the antenna load by: Vout = Ein ■ h = hEin ■ ua- 
The effective area and the magnitude of the effective height can be related through: \h\ ~ 

The polarization of h is aligned along the dipole axis 


y4e///(1207r) = Y'A^(Gam)/(4807r2) 
fiA ds h = hfiA, where \h\ is the magnitude of effective height. 

The full, complex transfer function T{uj) for the antenna, in principle, gives a complete 
description of the antenna (-|-cable) response and can be related to the effective height. 
This function T{uj) can be written as the product of the complex impedance of the (an- 
tenna-t-cable) Z{u}) multiplied by the complex height function h{uj), properly taking into 
account potential mismatches between the impedance of the antenna and the impedance of 
the attached cable: T = h{Zcable)/{Zantenna + Zcable), with 2'caWe=50ff+i0n. Both the magni¬ 
tude and phase of the effective height are determined directly from KUATR measurements. 
Zantenna IS determined from reflected power measurements on a HP8713C Network Analyzer 
(NWA). An independent check on the internal consistency of our Zantenna determinations is 
presented in Appendix 1. 

The magnitude of the effective height measured for a typical RICE dipole is given in 
Figure |^. The peak frequency in aii|^ (~600 MHz) is roughly consistent with the dimensions 
of the half-dipole (~15 cm); the effective height, as expected, also has magnitude of order 
10 cm. at the peak frequency. 

The phase variation of the effective height has also been measured as a function of 
frequency at KUATR and found to be negligible (put another way, we hnd that, to a very 
good approximation the experimental slope of the phase variation with frequency is linear: 
A0(0) 0.0310 (rad). Such a linear variation has no net effect on overall antenna response.). 

Given the magnitude of the effective height, the phase variation of the effective height, and 
therefore the magnitude and phase variation of the complex antenna impedance Z{ui), the 
complex transfer function T(a;) can now be calculated, and used to predict the expected 
waveform V'{t) observed in a RICE antenna in response to a transmitter signal. This is 
done by transforming the input impulse V(t) to frequency space, multiplying the function 
V(u) by the transfer function T(u) for both Tx and Rx (properly normalized), and then 
transforming back to the time domain to give V'(t). Figure [1^ shows the result of this 
exercise. Qualitatively, the after-pulsing observed in data is reproduced by our complex 
transfer function. We note that the actual expected signal shape observed for a neutrino 


®This peak frequency is shifted down in ice by n(a;)=1.78 at radio-frequencies. 
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Curve Fit (Fourier Sin/Cos Series Expansion) 



- Curve Fit - Real Data [ 

FIG. 12. Modulus of RICE dipole effective height (m.) as a function of frequency. Shown are 
data taken on the KU antenna testing range, with a Fourier fit overlaid. 

event requires knowledge of the details of the input signal, as discussed in the references 

00 - 

Dipole response has also been measured as a function of both azimuthal and polar angles. 
The polar angle response is observed to be well approximated by a cos^9 dependence (in 
power); the dipole response in azimuth is observed to be flat, as expected. 


C. Amplifier Gain Calibration 

Signals from the antennas are boosted by two stages of amplihcation, totaling between 
+88 dB and +96 dB of gain, depending on the channel. For a RICE waveform containing only 
(“unbiased”) thermal noise, the total power in a frequency bandwidth B can be calculated 
from the discrete Fourier transform of the waveform as Pnoise = kTB (we check several 50 
MHz-wide frequency bins from 250 to 500 MHz for this calculation). Since the total noise 
power in this band at the input to the antenna can also be written as a sum over the rms 
voltage measured in each frequency bin: P<y> = we can also write P<y> = PnoiseG, 

where G is the overall gain of the system. Thus, based on the rms voltage < H > of the 
8192 samples contained in these “unbiased” waveforms, the gain of the amplihers can be 
calculated in situ. The ampliher gain measured this way is flat up to the bandwidth limit 
of the oscilloscopes (500 MHz); direct laboratory measurements of the ampliher gain using 
a network analyzer are consistent with hat response up to 750 MHz. Most(>90%) of the 
amplihers are stable to 1-2 dB over the course of data-taking thus far analyzed. 
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FIG. 14. Comparison of expected (solid curve) vs. measured (histogram) Tx^Rx 
signal strength for one transmitter broadcasting to 16 receivers. Vertical scale is Re- 
turn-Power/Transmit-Power, in dB. No data is shown for the receiver (channel 2, top row) in 
the same hole as the transmitter being used for this test, due to possible cross-talk effects. Cor¬ 
rections for the measured roll-off of surface amplifier gain for some channels above 700 MHz have 
been only approximated. 


D. Full Circuit Amplitude Calibration 

In the final step of our amplitude calibration, the antenna response to a continuous 
wave (CW) signal broadcast from an under-ice transmitter is measured in situ. This test 
calibrates the combined effects of all cables, signal splitters, amplihers, etc. in the array. 
A 1 milliwatt (0 dBm) continuous wave signal is broadcast through the transmit port of 
an HP8713C NWA. The NWA scans through the frequency range 0—i>1000 MHz in 1000 
bins, producing a 0 dBm CW signal in each frequency bin. The signal is transmitted down 
through ~1000 feet of coaxial cable to one of the hve under-ice dipole transmitting antennas. 
The transmitters subsequently broadcast this signal to the under-ice receiver array, and the 
return signal power from each of the receivers (after amplihcation, passing upwards through 
receiver cable and fed back into the return port of the NWA) is then measured. Using 
laboratory measurements made at KUATR of: a) the effective height of the dipole antennas, 
as a function of frequency (previously described), b) the dipole Tx/Rx efficiency as a function 
of polar angle and azimuth, c) cable losses and dispersive effects (cables are observed to be 
non-dispersive for the lengths of cable, and over the frequency range used in this experiment), 
d) the gain of the two stages of amplihcation as determined from RICE data acquired in 
situ by normalizing to thermal noise Pnoise = kTB = < 1/^^^ > jZ, summing over all 

frequency bins, and e) hnally correcting for 1/r^ spherical spreading of the signal power, one 
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can model the receiver array and calculate the expected signal strength returning to the input 
port of the network analyzer. This can then be directly compared with actual measurement. 
Such a comparison, as a function of frequency, is shown in Figure for one transmitter 
(97Tx3). Below 200 MHz, the attenuating effect of the high-pass filter is evident. From the 
Figure, agreement is observed to be ~3-6 dB (in power) for the full circuit gain. Note that 
no correction for ice absorption has been made, given the small scale of the array.f^ 

Similarly, for each frequency bin we can calculate the difference between calculated vs. 
measured full-circuit gain. Figure [1^ shows the deviation between the calculated gain minus 
the measured gain, for several data runs. Included in the Figure are each of the 500 one MHz 
bins between 200 MHz and 700 MHz, for three transmitters. This Figure therefore shows 
the average deviation between model and measurement over that frequency range; for the 
five currently functional transmitters, the mean differences between the expected and the 
measured gainj^are —0.6 ±0.6, —0.8 ±0.6, —2.3 ±0.5, —3.4 ±0.6 and —2.8 ±0.6 dB. Within 
the “analysis” frequency band of our experiment (200 MHz - 500 MHz), our quoted level of 
uncertainty in the total receiver circuit power is ±6 dB; this value is commensurate with the 
width of the gain deviation distributions. On average, however, the calculated gains quoted 
above are well within these limits. Note also that, for all of the transmitters, the measured 
gain is actually higher than expected. To calculate eventual upper limits on the neutrino 
flux, however, we will use the more conservative calculated full-circuit gain. 


V. MONTE CARLO SIMULATION OF DETECTOR PERFORMANCE 

In order to check our understanding of the timing and gain calibrations of our experiment, 
we have written a Monte Carlo simulation of the RICE array. This simulation allows us to 
study the expected response of the radio receiver array to either a Cherenkov signal (as 
generated by a true z/giV —»• eN' charged-current event) or a random noise coincidence such 
as those expected to dominate our backgrounds. The simulation also checks our expected 
timing resolution (2 ns, as determined by measuring the channel-to-channel residuals for 
pulsed transmitter events in data) as well as study the effects of the overall gain uncertainty 
discussed previously (6 dB per channel, max, obtained by examining the residuals between 
calculated vs. measured gain for continuous wave signals from 200-700 MHz). Expected 
vertex (spatial) and angular (pointing) resolution, as well as neutrino energy resolution can, 
in principle, be assessed with the simulation. 


lONor have corrections been made for possible AMANDA cable “shadowing” in the same ice-hole, 
which is evidently not a significant effect. 

i^This is calculated as A(G'ca«c - G^eas) ± (^G,aic-Gm.as^ where /^{Gcaic - Gmeas) is the mean of 
each of the distributions shown in the Figure, and crc^^i^-Gmeas is the error in the mean, given by 
the r.m.s. of the distribution itself divided by the number of points in each distribution. 
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A. Monte Carlo Event Generation 


Neutrino interactions are generated uniformly in incident angle, with vertices uniformly 
populating the region: |x| <1 km, \y\ <1 km and 0 < z <1 km. After specifying a source 
location and incident neutrino energy, the simulation generates receiver times (smeared by 
our timing uncertainty of ~2 ns) and voltages at the input to the antenna, based on GEANT 

0 


simulations of radio signals in ice 


In addition to the voltage due to the neutrino 


interaction, thermal noise is also simulated at the input to the antenna. Antenna response, 
as a function of frequency, is modeled as described above; knowing the amplifier gain and 
cable losses, the voltage recorded at the surface is calculated, and smeared by the gain 
uncertainty of ±6 dB (±3 dB in voltage). With simulated times and simulated voltages, 
simulated events are then reconstructed with the same software used for data. 


B. Checks of the Simulation - MC vs. Data Tx Depth Reconstruction 

As a first check of the simulation, we have compared reconstructed transmitter source 
depths in data vs. Monte Carlo simulations. The source locations reconstructed from data 
collected when a pulser was connected to transmitter 97Tx3 were compared to 97Tx3 simu¬ 
lations, as shown in Figure 0- Vertices are reconstructed using the analytic 4-hit technique 
described previously. The simulation adequately reproduces both the location (to within 1 
m) as well as the width of the source distribution. 


C. Expected Vertex Depth Resolution Dependence on Depth 

As a further example of the utility of the simulation. Figure shows the expected vertex 
depth resolution (essential in discriminating surface sources from in-ice sources), as a function 
of the true source depth. Not unexpectedly, the resolution is best when the source vertex 
is close to the array and the event geometry is best resolved. Deep sources are increasingly 
difficult to pinpoint; additionally, the reconstruction software tends to reconstruct source 
vertices that are closer than the actual source {zreconstructed — Ztrue <0); i.e., vertices tend to 
be pulled closer to the array. This bias is under study. 


D. Angular Resolution 

Good angular resolution is essential in discriminating up-coming from down-going 
sources. Additionally, some physics analyses (e.g., coincidences with Gamma-Ray Bursts) 
are limited by the ability to point back to the recorded sky location of the GRB. For neu¬ 
trinos that interact very far from the array, the angular resolution is expected to be poor - 
the further the source point, the greater the inability to measure the Gherenkov wavefront, 
and the more the array “looks” like a single space-point relative to the interaction point. 
Figure |T^ displays the angular resolution for an ensemble of 10 PeV neutrinos interacting 
within 1 km in the ice below the array; we have required that at least four of the simulated 
receiver voltages exceed threshold and can therefore be used in vertex reconstruction. For 
approximately half the events, the angular resolution is about 10 degrees. The long tail in 
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FIG. 16. Comparison between reconstructed depth of transmitter (97Tx3), from pulser data 
taken at the South Pole, and Monte Carlo simulations of 97Tx3 pulser events. Times are smeared 
in the simulation by a Gaussian of width 2 ns. 
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A(z(reconstructed) - z(true)), meters 


FIG. 17. Deviation between reconstructed and true interaction depth {zAhit), 
action depth, from Monte Carlo simulations. The reconstructed vertex here was 
the analytic, 4-hit vertexing algorithm. The size of the squares corresponds to 
reconstructed events in the Monte Carlo simulation. 


vs. true inter- 
obtained using 
the number of 
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this distribution is due to distant source points which have correspondingly poorer pointing 
resolution. 



d0 (true-reconstructed direction) 


FIG. 18. Monte Carlo prediction for RICE array angular resolution, in units of degrees. Events 
{El, =10 PeV) have been simulated over the region within 1 km below and around the current 
array. 


E. Energy Resolution 

Having reconstructed a vertex location, the distance from each antenna to that vertex is 
determined. Having reconstructed the event geometry (i.e., the Cherenkov cone £t), we know 
the angular deviation of each antenna off the cone. The greater this angular deviation, the 
weaker the antenna signal strength. Finally, given the voltages recorded on each antenna, 
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we have enough information to make an estimate of the incident neutrino energy. For 
each antenna, the inferred value of Ey is calculated; we define Ey^reconstmcted as the simple 
average of the Ey estimates obtained from each antenna. Fig. p!9| shows the expected energy 
resolution {logio{Ey^reconstructed/Ey^generated)) for 10 PeV ueutriuos generated uniformly over 
distances within 1 km of the array. The resolution is obviously superior in cases where the 
hit antenna multiplicity is large {N^it >10) vs. cases where the hit multiplicity is small.[^ 


E(v) resolution: Nhit>10 {-) / Nhit<=10 {- -) 



FIG. 19. Histogram of number of events vs. logarithm (base 10) of Reconstructed/true neutrino 
energy, as obtained from neutrino interaction simulations. The sample is divided into high hit 
multiplicities (>10) and low hit multiplicities (<10). The energy resolution is obviously better for 
events with higher hit multiplicities. The simulated event sample is the same as that of the previous 
Figure. 


F. Conical vs. Spherical Source Event Geometries 

Perhaps most important in discriminating background from true neutrino-induced events 
is the characteristic conical Cherenkov distribution of energy in the latter case. Background 
(from the surface, e.g.) typically consists of spherical waves due to transient sources. We 
measure the consistency of a given event with either a conical or a spherical source by a 
“trial-and-error” procedure similar to the “grid”-based vertex-finding algorithm: given a 


^^In making this plot, the voltages on the antennas were smeared by the resolution on each channel 
(described previously); thermal noise has also been simulated. 
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reconstructed source location, we find the Cherenkov cone orientation most consistent (in 
terms of a minimum xLne variable) with the observed channel-to-channel voltages. Spherical 
sources should therefore correspond to large values of xLne since a conical source would 
typically have a smaller number of large-voltage signals than a spherical source. This is 
illustrated in Figure |^, where we compare the xlone distribution for: a) general, 4-hit 
RICE triggers, taken from data, b) calibration pulser events, taken from data, c) Monte 
Carlo simulations of a pulser source (97Tx3) emitting spherical waves, and d) Monte Carlo 
simulations of neutrino events producing Cherenkov cones. The separation between conical 
vs. spherical source geometries is evident from the Figure, as well as the rough agreement 
between data and the simulation of a spherical source near the array. 


Comparison of data vs. MC cone-fit chisquares 



FIG. 20. logio{x^) of fh of recorded voltages (for both data and simulation) for spherical vs. 
conical geometry sources. The separation between conical and spherical sources, as well as the 
agreement between simulation and data for spherical sources, is evident from the Figure. 


VI. SUMMARY 

Basic calibration of both the time and amplitude response of RICE radio receivers has 
been made, relying primarily on data taken in situ. A Monte Carlo simulation has been 
written which reproduces the gross features of those calibration data. The calibration of the 
detector is sufficient to allow limits to be placed on the incident high-energy neutrino flux. 
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APPENDIX 1: Check of Antenna Impedance 
We can check the internal consistency of our Zantenna measurements by constructing an 
appropriate Green function to give the antenna response to a given input. Impedance is, in 
effect, a Green function - namely, an analytic function in the complex uj plane, reflecting the 
constraint of causality in the time domain. 

The Zeroes and Poles Expansion of the impedance of a network implements the essential 
analytic features. It states: The impedance of a network containing any number of L, R, C 
elements, including antennas, is a rational function of ut, which can always be factored into 
simple zeroes and poles ||22|| : 


Z{u;) = , 


(a; - a;o)(a; - a;i)(a; - a;2).. 

iuj(uj — uJo)(uj — a;()(a; — taf}... 


( 1 ) 


The expansion follows from the dehnition Z{uj) = V{uj)/I{uj) and the expansion of V{u!) 
and I{uj) into polynomial ratios, which can be factored. The pole at a; = 0 is just the static 
capacitance of the device. The poles and zeroes are symmetric about the imaginary axis and 
all lie in the lower ut half-plane, consistent with causality. 

Since all the poles and zeroes are simple and isolated, an equivalent ansatz is a sum of 
familiar “Breit-Wigner” resonances. 


iZH = E 

J 


—iajujujj 
a;2 — cuj J- iVjui 


( 2 ) 


Here ujj is the value of a resonant antenna frequency, and Fj is the width; the constant aj 
measures the strength. 
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Check of impedance ansatz with data 


A Hewlett-Packard network analyzer HP8712C measured real and imaginary reflection 
coefficients P^, P* for antennas in the lab. Tests were made after calibrating out cable, 
connector and lead effects. The reflection coefficients are related to the complex impedance 
Z{uj) by 


Re{Z{uj)) = 


(i-p,)2 + r2 


and Im{Z{ijj)) = 


2P. 


(i-n)2 + rr 


( 3 ) 


Variations from antenna to antenna were observed to be small, typically 10% or less and 
considered smaller than other variables in the calibration chain. 


Figure ^ shows the real and imaginary parts of the antenna impedance measurements 
over the approximate range 50MHz — 5GHz. This range substantially exceeds the range 
over which the antenna impedance needs to be known, due to cutoffs at the low and high 
frequency ends from the high-pass filter and cable losses, respectively. These effects are 
well-characterized and discussed elsewhere. The Figure shows real and imaginary parts as 
“bumps” and wiggle-excursions consistent with Breit-Wigner physics. 

We tested our understanding of the antennas by fitting the real part (only) of the 
impedance to a sum of Breit Wigner functions. This is straightforward because the res¬ 
onant frequency is located near the bump maximum and the width can be estimated from 
the graphs. Fitting one-half of the data is the same procedure as fitting a dispersion relation, 
namely 


iZ{uj) 


1 

27ri 


/•+°° , Im{iZ{u}')) 

/ du - - — 

J—oo CO — CO' 


Compared to standard dispersion relations, real and imaginary conventions are reversed in 
impedance, with the absorptive part of impedance (resonant bumps) defined as real-valued. 

Once the real part was fit, it was used to predict the imaginary part of the impedance 
(Fig. ^ gray lines). The predictions agree with the data and serve as a check of the entire 
procedure. 


30 







FIG. 21. Real (positive only) and imaginary (both positive and negative values) antenna 
impedance data obtained from network analyzer measurements (dots), as a function of frequency 
(MHz). The vertical scale is linear and in absolute units of 50 ohms. After performing a multiple 
Breit-Wigner fit to the real part of the impedance function Re{Z[ijo)), the imaginary part of the 
impedance is predicted directly from the real part by the dispersion relation. Agreement between 
the calculated and the measured values of Im{Z) is satisfactory. 
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